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Abstract 

We calculate the inclusive decay width of T — > l + l~ . Then we get the ratio R Tfl = T[T — > 
r + T~]/r[T —7- to 0(a) and 0{a 2 s ) within the Standard Model(SM). Comparing with the 

recent Babar's data R Tfl = 1.005 ± 0.013 ± 0.022, we hnd that SM prediction R Tfl is not consistent 
with the experimental data in the error bar. The discrepancy is about 1.25<r. So leptonic decay 
of T may be a possible signature of New Physics(NP). We present a better approach to test the 
Standard Model, R Tfl {E soft ) = T[T -)• t+t" + X]/T[T fj, + fi~ + X]\ Ex<Eaoft . After resumming 
the large logarithms , we get R T ^{E so ft) at the precision level of 0.1%. It can be compared with 
experimental data more precise. We also consider the impact of R T ^(E so ft) and R T ^ from light 
Higgs h and pseudoscalar Higgs Aq. 
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I. INTRODUCTION 



Although the Standard Model(SM) of particle physics describes the interactions of ele- 
mentary particles very successfully, it is believed that SM is not the final theory and there 
should be New Physics (NP) beyond SM. So the hunting of NP is one of the hottest topics 
for theorist and experimentalist. The B factories gave a very clear channel to test SM, just 
as T(35) — > T(1S)tv + tv~ , T — > l + l~ (I = r, //). Recent Babar measured the ratio 

elq 

Br\T -> t + t~] 

R ^ = 1TF? — ^ = 1.005 ±0.013 ±0.022, 1 

Br[T — > /U+^^J 

where branch ratio Br[T — > t + t~ (fi + fi~)\ is corresponded to inclusive decay width. The 
final states radiations(FSR) effects due to photon(s) and gluon(s) are taken into account in 

nn 

MC generator. The Leading Order (LO) SM prediction of R Tfl is 0.992 |3|, |4f. It is consistent 
with experimental date. Then Babar claimed "No significant deviation of the ratio R Tfl from 
the SM expectation is observed" . 

Theoretically, the high order corrections of the ratio R Tfl should be taken into account. 
The SM predictions should be compared with experimental data beyond tree level. At the 
same time, R Tll is sensitively on the coupling of h(Ao)bb and h(Ao)l + l~ within NP. It is 
an excellent probe for the new Higgs interactions in some NP Model, where the coupling 
of Higgs bb and Higgs is enhanced Then we should calculate the ratio R Tfl and 
compare with the experimental data to test SM or hunt NP. 

There are some theoretical and experimental works related with leptonic decay of T. The 
Quantum Chromodynamics(QCD) corrections of T — > have been calculated to two- 
loop [6]. We have calculated T decay to charm jet [7]. The leptonic decay of vector bosons 
has been calculated to Next-to-Leading Order(NLO) in Quantum Electodynamics(QED) 



,8 



The CLEO got the ratio R Tfl = 1.02±0.02±0.05 in 2006 |9j. The MC simulation of T -> l+l 



has been studied, where large logarithms 
A is also introduced in those processes 



rave been resummed|l0l|. The pseudoscalar Higgs 
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12j . Babar has searched for a light Higgs 



boson A in the radiative decay of T(nS) — > 7v4 , A — > l + l for n — 1, 2, 3. They found no 
evidence for such processes in the mass range Q.212GeV < Mao — 9.3Gel^ and no narrow 



structure with 4.03Ge^ < M T+T - < 10.10GW [l3 



Q. 



In this paper, we calculate the inclusive decay width of T — > l + l . Then we get the precise 
prediction within SM. We also consider the impact from light Higgs h and pseudoscalar Higgs 
A . 
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FIG. 1. Part of the Feynman diagrams of T — > l + l within Standard Model. 



II. SM PREDICTION 



The LO QED Feynman diagrams of T — > l + l are shown in Fig. [TJ Followed the process 
of T cc in Ref. [3], we can get the LO amplitude and decay width of T — > l + l~ , 



M LO [T^l + r] = x l^a\R(0)\l/l, 



Tlo[T 4ii| = ^ , (2) 

where r; = Mf/My, \R{0)\ is the radial wave function of T at origin, e is the polarization 
vector of T. If expanded with ri, we can get 

r LO [T -+ l+n = 4| ^[ a2 (1 - 6r? + O (rf)) . (3) 

We take into account the NLO QED correction here. The renormalization of lepton and 
b quark wave function, and electron charge should appear. We use D = 4 — 2e space-time 
dimension to regularize the divergence. On-mass-shell (OS) scheme is selected for Z 2 ui) and 
modified minimal-subtraction ( MS ) scheme for Z P : 



"®=--k 



1 2 , 47T/X 2 





euv em M z 



^ e MS = f (3 + ^)f— -7* + lrW), (4) 
bit 3 Ve uv / 

where fi is the renormalization scale, 7^ is the Euler's constant, / = b,l, and Qf is the charge 

of fermion / in unit of electron charge. The factor 3 + y is from the charge and color factor 

of three flavor lepton e, fi, r (3 x 1) and four flavor quark u, d, s, c (2 x 3 x (1/9 + 4/9)). If 

we ignore the self energy of photon and the renormalization of a, the NLO QED correction 

is just replaced 4a s /3 with a from T — > ccQ]. 



TABLE I. The numerical decay width of inclusive processes T — > l + l (I = r, /x) in unit of jnTgy^ 
and R TLL within SM. 





t[ T ] m 


Rt/j, 


LO 


2.822 2.844 


0.992 


NLO QED 


2.777 2.798 


0.993 


NLO QED, l + l-gg 


2.780 2.836 


0.980 


NLO QED&QCD, l+l~gg 


1.743 1.791 0.973 ±0.001 


Babar 


- 1.005 ± 0.026 



For the corrections due to gluons in the final state are considered in experimental Monte 
Carlo, we should consider the QCD processes T — > l + l~ + gg. We also consider NLO QCD 
corrections to the decay width of T — > l + l~ , which give a factor of 1 — 4:Cfa s /ir to suppress 
the LO decay width js]. 

In numerical calculation, the parameters are selected as: 

M e = 0.5110MeV, M d = O.OOMeV, M u = O.OOMeV, 

M M = 0.1057GeV, M s = O.lOGeV, M c = 1.30GeV, 

M T = 1.7768GeV, M b = 4.73GW, a = 1/132.33. (5) 

Here = My/2. The renormalization scale \x is selected as fi — Mx, and the fine structure 



constant a is calculated with the program alphaQED.f [15(].The numerical r[r(u) and R T/X 
are listed in TableJU The LO prediction of R T ^ is 0.992. It is used in Ref. |l|, 12], where 
claimed "No significant deviation of the ratio R Tfl from the SM expectation is observed". 
But the QCD corrections should suppress the SM prediction and drive R Tfl away from the 
experiment data. 

We should calculate the uncertainty for the theoretical prediction. As an order estimate, 
one can get R^® ~ 0((a/ir) ). For the NLO QED corrections have been taken into account, 
the uncertainty from higher order QED contributions is 0(a 2 /ir 2 ) ~ 6 x 10~ 6 . In the same 
way as QED, the uncertainty from higher order QCD contributions is 0(a? s /'K 2 ') ~ 3 x 10~ 4 . 
Z can contribute to T — > at tree level through replacing photon with Z. We can get 

M z LO [T -» /+/-] = r[(4sin 2 ^-l) ^7 5 H ffil 
Ml [T^l+l-} Iz I 41 ' U 



and 

= M 2 (3-4sin 2 ^) 

h 16 (M 2 - M§) (1 - sin 2 O w ) sin 2 % ' U 

Here f z ~ —M\jM\ ~ — 10~ 2 . The vector current term should change the LO amplitude 

in Eqf2]by a factor f z (— 1 + 4 sin 2 9w) ~ 10~ 3 . It is just like replacing a with a(l — f z (l — 

4 sin 2 6w)) in Eqj2], but it is a global factor for three lepton at LO. Then the uncertainty from 

Z of R TI1 should be 0(f z (l - 4 sin 2 9 W ) (R^ D - R™)) ~ C(10- 6 ). Here superscript Q^D 

means NLO QED has been taken into account. The axial vector current is not coherent with 

the vector current in Eqf2J It change the width with a factor O(M^fM^) ~ (9(10~ 4 ) and the 

ratio with a factor 0[M\M^ / 'Mf) ~ 0(1O~ 5 ) only. Compared with the Z contribution at 

tree level, the other contributions from weak bosons should be suppressed by a/n or more. 

Then the weak contributions from W ± ,Z,H can be ignored safely. Within SM, it should 



be considered that T — > 7775, where r\ h — > l + l is followed 12 1. The energy of 7 is about 
70 MeV in T ^ 7 r/ b and Br[r] b Z + r(+7 so/t )] ~ 10" 8 Q 



17j . For T — > 777^ is a P wave 



process, we can estimate Br[T — > 777^] through 



r [T -> 7 77 6 ] M w (M r - M, 



r [ J/V -> 7 77 C ] V e c 7 V m t (M J/V , — M Vi N 1 
Then 5r[T 7%] ~ 10~ 5 . So 73r[T -> 7775] x Sr^ ->• / + /"(+7 so /*)] ~ !0~ 12 - Thi s can be 
ignored safely. The uncertainties of R Tfl within SM are listed in TabJLTl Then SM prediction 
is 

R T ^ = 0.973 ± 0.001. (9) 

Compared with Eq.([T]), it is not consistent with the experimental data in the error bar. The 
discrepancy is about 1.25a. 

The QCD contributions have been taken into account in EqQ. It is difficult to 
measure. So we present a better approach to test the Standard Model, R Tll (E so f t ) = 
T[T t+t- + X]/T[T + X\\ Mx<E If we select E soft ~ 5GeV, T[T l+l~ + 

gg\\M x <E soft is less than T[T — > / + Z~]/1000, then the impact on R T ^(E so f t ) is less than 
2 x 10~ 5 , but the large logarithms appear 

, 4£ 2 4M 2 , , 

^tA^^ft- (10) 



M{ Mf 



We resum the large logarithms with YFS resummation scheme 10, ll8|] 



—a ( , , 2E 1 , lnr; 7r 2 , 

r = _N (lnri + 1)ln _^ + _i__ + 1 |. 



TABLE II. The uncertainties of R TI . within SM. 







Order 


Numerical 


QED 
QCD 

Z{W ± ,H) 
Vb 


a 2 /vr 2 
a 3 s /TT 3 

M\MflM% or aM 2 /(Mf vr) 
Br[T -> 7%] x £r[% -> Z+Z - 


6 x 10~ 6 

3 x 10~ 4 

4 x 10~ 6 
1 x 10" 12 


Total 


0.001 


r>SM 




1 


0.973 ± 0.001 



The resumed results are 

r% = e Y r LO , 

FnIo = (e Y -l-Y) T LO + T QED . (12) 

After the large logarithms are resummed, we get R Tll (E so ft) with a soft cut at the precision 
level of 0.1%. The numerical T T ^(E so f t ) in unit of -^g^ and R T ^(E so f t ) with different 
energy cut E s are listed in Table JLTTl The dependence of R T ^(E so f t ) on the soft cut E s is 
shown in FigfSJ If we select E s = 0.2GeV. Including the uncertainty, the ratio is 

R T ^(0.2GeV) = 1.0628 ± 0.0011. (13) 

The effect of QCD is very weak in this channel. R Tfl (E so f t ) can be compared with experi- 
mental data more precise. 

III. IMPACT FROM NP 

NP may play a role in the discrepancy between theoretical prediction and experimen- 
tal data of R TfJi in Eq.flH]) and Eq.((T]). We only consider the scheme of light Higgs h and 
pseudoscalar Higgs A here. 

The Feynman rules are ieMfCh/ (2Mw sin 9w)ff and — eMfCAo/{2Mw sin 9w)fj 5 f for 
hff and A ff vertex respectively, here f — l,b. Cao(Ii) are different in the special model, we 
consider them as parameters. For it is IR finite which A (h) involved in T — > 7 so /t/ + /~, so its 
contributions are suppressed by E s /M b ~ 4 x 10~ 2 when compared with virtual processes. So 
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FIG. 2. The dependence of R TfJi (E so ft) on the soft cut E s within SM. 

I /?foil 2 

TABLE III. The numerical decay width of processes T — > = r, fj,) in unit of jgr^gr^ and 

R Til {E so f t ) within SM. E s = 0.1 means the soft cut is O.lGeV . 





r[r] m 


RrfiiEsoft) 






LO 


2.8221 2.8444 


0.9922 


LOYFS 


E s 


=0.10 


2.7277 2.4925 


1.0944 


NLO 


E s 


=0.05 


2.6744 2.3932 


1.1174 


NLOYFS 


E s 


=0.05 


2.6768 2.4272 


1.1028 


NLO 


E s 


=0.10 


2.6954 2.4678 


1.0922 


NLOYFS 


E s 


=0.10 


2.6970 2.4916 


1.0824 


NLO 


E s 


=0.20 


2.7158 2.5411 


1.0688 


NLOYFS 


E s 


=0.20 


2.7168 2.5564 


1.0628 


NLO 


E s 


=0.45 


2.7385 2.6236 


1.0438 


NLOYFS 


E s 


=0.45 


2.7389 2.6312 


1.0409 



we ignored the real processes and included the virtual processes only when we considered the 
impact of A (h) to R TfJi (E so f t ). The Feynman diagrams are shown in FigJ31 The Feynman 
diagrams which exchange A (h) between bb are ignored for it should not change the ratio 
R Tfl . Compared with T[T — > t + t~], the impact of A (h) to T[T — > is suppressed by 

M^/M^, for the A (h)l + l~ coupling and the spin flip between l^ v l~ and are both 
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4 NP Box 2 NP Vertex 



FIG. 3. Part of the Feynman diagrams of T — > l + l which Aq(K) involved. The Feynman diagrams 
which exchange Aq{K) between bb are ignored for it should not change the ratio R TfjL . 
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FIG. 4. The Aq{K) impact on T — > t + t~ as a function of Mj^qqa. The Aq(K) impact on real 
contributions ignored for it is suppressed by E s /Mb and T — > is ignored for it is suppressed 

by M^/M^. The Feynman diagrams which exchange AqQi) between bb are ignored for it should 



not change the ratio R T , 



proportional to M\. So we ignore the contributions for T — > . Only T — > A^* — > t + t~ 
and T — > h*j* — > t + t~ are taken into account. The contributions with AqAq, h*h*, or A* Q h* 
in the loop are zero for J PC . The numerical result of the Ao(h) impact as a function of 
M^offt) from the loop Feynman diagram is shown in FigJH The Aq(K) impact on R Tfl is 
RL°T m W[t]/T lo [t\. 

If we consider the R T ^, we should include the real correction too. If we select 10.3GeV < 

M Am < 10.6GeV, T ao [t}/T lo [t] 4 x l(T 6 Ci + 5 x lCr 10 C%, and T h [T\/Y LO [r\ ~ 

3 x lCT 6 ^ + 8 x 10- 10 Ct 
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IV. SUMMARY 



In summary, we calculate the inclusive decay width of T — > l + l (I = r, yu). then we get 
the ratio R Tfl = T[T — > r + r~]/r[T to 0(a) and 0(o£) within SM. Compared 

with the recent Babar's data R TfJL = 1.005 ±0.013 ±0.022, we find that SM prediction R Tfl = 
0.973 ±0.001 is not consistent with the experimental data. The discrepancy is about 1.25a. 
So leptonic decay of T may be a possible signature of NP. We present a better approach 
to test the Standard Model, R TIM (E soft ) = T[T -> t+t~ + X]/T[T -> + X]\ Ex<E . 

After resumming the large logarithms, we get R Tll (E so f t ) with a soft cut at the precision 
level of 0.1%. The effect of QCD is very weak in this channel. It can be compared with 
experimental data more precise. We also consider the possible solution, light Higgs h and 
pseudo scalar Higgs A . To clarify the discrepancy, more work should be done by theorist 
and experimentalist. 
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